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Abstract 



arating gluonic fields in background and quantum fields and neglecting an- 
harmonic contributions. Mesonic states with nonperturbative gluonic compo- 



X 

nents are constructed. These states are invariant under gauge changes of the 
background fields and form an orthogonal basis. The effective hamiltonian 
is diagonalized in this basis in a systematic 1/m- and short distance expan- 
sion. For very heavy quarkonia, we obtain an effective potential similar to 
the phenomenological funnel potential. We compare our method to 2 nd order 
perturbation theory in the background fields and demonstrate its applicabil- 
ity even for the relatively light charmonium system. The results to order 1/m 
for pseudoscalar meson masses and wave functions are shown and compared 
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with those of the Cornell model. 
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I. INTRODUCTION 



In phenomenological hadron models nonperturbative gluonic effects are accounted for 
in a variety of ways. In non-relativistic potential quark models Q, an effective interquark 
potential is assumed to result from them and the effective hamiltonian is diagonalized in the 
Hilbert space of quarks only In flux tube models one does not completely eliminate the 
gluons. Their net effect is to generate a color flux tube between quark and antiquark that 
binds them. In bag like models || nonperturbative gluonic effects appear in the guise of the 
bag constant. They are taken to generate a vacuum pressure which counterbalances that of 
the (perturbative) quarks in the interior of the bag. Gluons in the bag are assumed to be 
perturbative and responsible for a hyperfine interaction. All these models are surprisingly 
successful in describing the hadronic spectrum. Except for a few exotic states there seems 
to be no need for hybrid states, bound states of quarks and gluons. Nevertheless, even one 
gluon exchange, leads to intermediate states, where the (constituent) quarks are not in a 
color singlet representation and one may wonder, whether the nonperturbative part of their 
interaction (for our purposes everything that is not one gluon exchange) is really as color 
blind as is generally assumed. A better understanding of the relation among the various 
model parameters and the structure of the nonperturbative ground state would also clearly 
be of interest, and there have been attempts to find a relation between the bag pressure 
and gluon condensation, and the confining potential in nonrelativistic models and the string 
tension obtained in lattice calculations f|. 

In this paper we study the influence of a nontrivial gluonic ground state on the structure 
of heavy quarkonia. The nearly nonrelativistic nature of these mesons makes them ideal 
probes of ground state properties. Their large mass allows for a systematic expansion of the 
interaction at small distances and we gamble that the lowest dimensional (gluonic) conden- 
sate suffices for a rudimentary description of the nonperturbative vacuum structure in this 
case. We estimate the matrix elements of the hamiltonian in an approach very similar to the 
one used in QCD-sum rules [[| using background fields to describe nonperturbative gluons 
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and parametrizing their vacuum matrix elements. In contrast to the approach taken by 
Voloshin |J and Leutwyler [[!]] , we do not treat the nonperturbative part of the hamiltonian 
so obtained as small compared to the (perturbative) coulomb interaction. This is probably 
only the case for quarkonia heavier than botonium 0. It has been previously suggested that 
the inclusion of the finite correlation of the gauge invariant correlator 

SjsSUO*) =< I Tr{F, p (x)S(x, 0)^(0)5%, 0)} | > (1) 

will strongly diminish the importance of couplings to the background fields, allowing again 
for a perturbative treatment of them [|]|§. In eq.fll]) 

S(x,0) = Pexp(z f 1 dtx^AJxt)), 
Jo 

A, = gT a A^ 

S is the color transport operator required for gauge invariance; T a are the generators of 
the gauge group in the fundamental representation. The correlation length of should 
be no larger than that provided by the mass of the lowest physical state which contributes. 
In a purely ggluonic theory this can at best be that corresponding to the lightest glueball. 
According to lattice estimates, the lowest glueball mass is probably larger than 1500 MeV 
||10|| , and the correlation length of ([J) therefore is probably less than 0.13 fm. The value of 



0.22 fm extracted by Di Giacomo et al. JTTj] from a simulation on a modestly sized lattice 
is probably compatible with this estimate, once errors and finite size effects are taken into 
account. If one accepts the arguments used for instance in ref. the nonperturbative 
effects due to condensates in heavy quarkonia would be exceedingly small and of little 
importance. To reproduce even the first low lying states of the quarkonia spectrum additional 
nonperturbative effects, such as a confining potential are required in the approach of ||. 

We wish to point out however, that the correlation function ([]]) is not the relevant one 
for heavy quarkonia, because it describes the propagation of gluonic modes in the absence of 
heavy quarks. This correlator has very little to do with the low energy interaction in heavy 



quarkonia: the heavy qq pair is only rarely destroyed to produce gluonium intermediate 
states - this is a Zweig rule suppressed process, necessitating relatively large momentum and 
thus having a relatively short correlation length. In heavy quarkonia there almost always is 
the qq pair around to screen the color charge of soft effective gluons being exchanged. We 
propose that the relevant correlation function is that of a composite color singlet (scalar) 
operator such as ^>F^ in a quarkonium state | M >, written schematically as 

g M (x -y)=<M\ (*F*)(a;)(*i5 , *)(y) | M > . (2) 

Essentially two different kinds of physical intermediate states contribute to this correla- 
tion function: 

1. Glueballs: in this cases the correlation functions factorizes into a decay matrix element 
of the quarkonium and the correlation function (|TJ) - the correlation length of this 
contribution is essentially given by the lowest glueball mass. It is the Zweig rule 
suppressed process we mentioned, and we will neglect this contribution for the reasons 
given above. 

2. Quarkonium states: they only contribute, if the quarkonia have admixtures, where the 
quark- ant iquark pair is in an octet configuration (with additional gluons forming an 
overall color-singlet meson). If only quark- ant iquark (color-) singlet configurations are 
considered, the contribution from the excited quarkonia to (0) would vanish by color 
invariance. Since the mass difference of the lowest quarkonium to the first excited 
excitation is of the order of 600 MeV, one estimates that the correlation function @ 
will have a correlation length of about 1/3 fm - much larger than one would expect from 
glueball intermediate states. As explained in section |T| we do approximately include 
such correlations by extending the basis Fock-basis for the description of quarkonia by 
states in which the heavy QQ-pa.ii is in an color-octet configuration, with an additional 
(soft) color-octet effective gluon to form the singlet. 

We calculate the short distance expansion of the hamiltonian matrix elements in a gauge 
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invariant basis of color singlet states ||12|| , including only terms up to order r 2 and 1/m. We 
then obtain the eigenstates by (numerical) diagonalization of this hamiltonian matrix and 
see how far we can go in this approach. The previous discussion implies that a realistic 
description of heavy quarkonia within this short distance expansion can only be obtained 
if the (background) gauge invariant basis of ref. is extended to include states in which 



the heavy quark-antiquark pair are in a color octet representation but coupled to a nonper- 
turbative gluonic background field ||[13| to form an overall color singlet. Diagonalization 
of the hamiltonian matrix leads to a "color-octet" component in the wavefunction of heavy 
quarkonia that takes account of the possibility of color exchange between the valence quarks 
and the background field. This component becomes rather large for higher excited states 
and its coupling to the "singlet" is the main reason for the distortion of the coulombic spec- 
trum in this model. In the pseudoscalar meson channel only the few octet states constructed 
in section [II] couple and a numerical diagonalization of the resulting system of differential 
equations is still quite feasible. 

Recently, calculations of heavy quarkonium annihilation rates have acquired a rigourous 



theoretical treatment through general factorization formulae obtained in the context of 



the effective field theory, NRQCD ]T5| . In the factorization formulae short- distance coeffi- 



cients, calculated in perturbation theory, are combined with long-distance matrix elements 
that can so far only be computed in lattice simulations. This approach however provided 
the solution of the long-standing problem of infrared divergences in he and Xci decays into 
light hadrons. The probability for the qq pair to be in a color-octet s-wave at the origin is 
an essential input in this approach. Since accurate lattice calculations of this probability are 
not yet avaiable, it is determined phenomenologically fl6 |. In the approach we propose, this 



probability can in principle be computed, since such octet configurations are included in the 
description of the meson from the outset. In this sense, our approach could be considered 
an inexpensive alternative to lattice simulations. 

In section |ITJ we derive the effective hamiltonian appropriate for the description of heavy 
quarkonia. We first obtain the effective Lagrangian to order 1/m by a Foldy-Wouthuysen 



17 1 transformation. Using the background field formalism and neglecting anharmonic 



quantum fluctuations one finally arrives at an effective Lagrangian |^9|,^0j that includes 
background fields up to 2 nd order and which, in the instantaneous approximation, gives rise 
to an effective hamiltonian that is accurate to order 1/m and r 2 . Since retardation effects are 
of order 1/m 2 , their consistent inclusion would require a much more elaborate treatment, 
which would only obscure the basic nonperturbative gluonic effects we want to elucidate 
here. 

Numerical diagonalization of this effective hamiltonian in the extended basis then yields 
(pseudoscalar) meson masses and wave functions. Our approach can in principle not be 
described by an effective potential in the singlet channel |],[n| because the elimination of 
the relative octet states would make it energy dependent. One can however obtain an energy 
independent effective potential for infinitely heavy quarks, where all terms of order 1/m (also 
the kinetic energy) can be neglected. This potential should be closely related to the static 
quark potential one extracts from the Wilson loop [^TJ. This potential, which we derive 
and discuss in section [TV], also gives us an idea how far the short distance expansion can be 
trusted. We show that although it is apparently linearily rising at intermediate distances 
Afm < r < .7fm, the potential is also compatible with an effective exchange of the Gribov 



type [|22 



In section [V] we compare our method to methods relying on 2 nd -order perturbation theory 
in the background field and to the phenomenological Cornell |31| potential. We first present 
results where all potential matrix elements of order 1/m are neglected. This greatly simplifies 
the calculations because only a few basis states couple, but pseudoscalar and vector mesons 
are degenerate at this level. We also give the results of a more complete calculation of the 
pseudoscalar quarkonia which includes the 1/m corrections. 

Section |VT| is a summary and discussion of our results. 
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II. GAUGE INVARIANCE AND BASIS STATES 



It is usually assumed that physical states are color singlets. For heavy mesons, where the 
non-relativistic approximation is adequate, one can represent quark and antiquark fields by 
2-component spinors. Gauge invariance of the state requires that color is parallel transported 
from the the quark to the anti-quark along some path. With a straight path the quark and 



simplicity we will restrict ourselves in the following to pseudoscalar mesons. The simplest 
gauge invariant pseudoscalar state is of the form 



where w a (x2) creates a quark with color a and spin a at f 2 - v does the same for an antiquark. 
We will refer to such gauge invariant states as singlet states since the quark anti-quark pair 
at vanishing separation is in a color singlet representation. 
We choose for the color transport operator 



the path ordered exponential (denoted by Pexp) of gluon operators A^(x) along a straight 
line from x\ to x%. It is then relatively straightforward to show that canonical anti- 
commutation relations for the quark and anti-quark operators imply that the singlet meson 
states (||) are orthogonal fl2|| . 

We can also construct gauge invariant basis states where the valence quark anti-quark 
pair is in a color octet representation at vanishing separation by coupling them to chromo- 
electric or -magnetic fields. We shall call these states octet states for obvious reasons. Since 
the chromomagnetic field B = B a X a /2 transforms as a pseudovector and the chromoelectric 
field E = E a X a /2 as a vector under rotations and according to the adjoint representation of 
the gauge group, we extend the basis for pseudoscalar mesons by the states 



anti-quark anti-commutation relations ensure that the basis states are orthogonal |fL2|j . For 




(3) 




(4) 





S 



I 21 )ei = J2 ^j-u ] a (x 2 )E(x 2 , t) ■ (x 2 - fi)T(f 2 , xi)v a (xi) I ft) , (6) 

and 

I 21) £2 = Y] rr ui(x 2 )E(x 2 , t) ■ (a a/3 x (x 2 - fi))T(f 2 , £i)^(£i) | ft) , (7) 

a/3 V27T0 

the summation over color indices being implied. 

The above states are seen to be mutually orthogonal and normalized by the canonical 
anti-commutation relations of the quark and anti-quark operators if we assume expectation 
values 

a 2 a 2 
(ft I ^—B ia B> b I ft) = -(ft I -f-E ia E jb I ft) 
47r 1 ' N 1 47r 2 1 ' 

= — W 6 (ft I -F» UC F' I ft) = —6 ij 6 ab 6 2 , (8) 
96 x 1 7T ^ 1 7 96 v ; 

and < i? >=< -B >=< EB >= 0, which are a consequence of the Lorentz- and parity- 
invariance of the vacuum state | ft). Its nontrivial nature is reflected in a non- vanishing value 
for (p 2 , which from QCD sumrule estimates should be close to (SSOMeV") 4 ||. If we neglect 
expectation values of all higher dimensional gluonic operators the states (0), (§),(]§) and (|7|) 
form a complete orthogonal basis for the valence quarks of heavy pseudoscalar mesons, while 
the singlet states (§) are only complete in this sense if one neglects the non-trivial vacuum 
structure altogether. 

We thus effectively truncate the Fock space by taking basis states whose gluonic sector 
contains at most 1 soft gluonic mode, i.e. we replace the low momentum gluonic modes 



with octet quantum numbers by zero-modes. We show in section |T| that such a severely 
truncated basis is adequate fo the problem at hand. 

Note also that Lorentz invariance of the ground state |ft > forces one to assume that the 
chromoelectric background field is either antihermitian or that the "E-states" have negative 
norm (see (§)). Both possibilities lead to a nonhermitian hamiltonian matrix, whose eigen- 
values in general are not real. Already the construction of basis states for the mesons on a 
non-trivial ground state indicates that one can at best hope to find a few stable mesons in 



this approach. We assume that a pseudoscalar quarkonium r\ is well described by a linear 
combination of the above basis states 



V)= E / i > A f(2,l)|2,l) M . (9) 

M=S,El,E2,B' 
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III. EFFECTIVE HAMILTONIAN 

We would like to use the basis constructed in section 2 to approximately diagonalize 
the QCD hamiltonian for heavy quarkonia. Since low lying heavy quarkonia are generally 
believed to be of relatively small size, the interaction should be dominated by the pertur- 
bative gluon exchange, which leads to a coulomb-like effective potential. The experimental 
spectrum deviates however noticeably from a purely coulombic one and can be reasonably 
well reproduced by the combination of a coulomb- and a linear confining- force. We wish 
to emphasize here that the truly long range part of the confining force (r > lfm) is not 



really tested by the observed quarkonium states |23|. Our conjecture is, that what has to 
be included in a systematic approach are the first short distance corrections to the (pertur- 
bative) coulomb force due to the non-trivial structure of the gluonic ground state. This was 
proposed previously but without any tangible results, because the corrections were found 
to be exceedingly large within the framework of 2 nd -order perturbation theory except for 
quarkonia beyond bottonium |],[7| . This apparent failure of an idea which we believe is quite 
well founded, has led us to reexamine the basic procedure used in the evaluation of these 
effects. 

In this section we outline the derivation of the effective hamiltonian for heavy quarkonia, 
whose matrix elements are correct up to order a,l/m and r 2 . To this order in the short 
distance and heavy quark expansion the nonperturbative aspects of the ground state can be 
described by the gluon condensate <fi 2 . 

The nonrelativistic approximation for heavy quarks is conveniently obtained from the 
QCD lagrangian 
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c(x) = -\f; v f^ + + 9T a \f a )i> - , (io) 



by a Foldy-Wouthuysen transformation W%. In terms of transformed quark fields 



ip — > exp(ry • D/2m)ip , 

—> ip exp(—rf ■ D /2m) , (11) 



the lagrangian is 



£nrqcd {%) = -\f; v F>? + ^7° A) - m)V 



+V>(^[7 " A A)] + L2 ^))^ + 0(l/m 2 ) 
= -~^Vr + #7° A, " ™)V + (12) 
-^^--^ + ^^— ^ + 0(l/m), 

where D = d — igV and E = 2) does not couple upper and lower spinor components. 
They are only coupled in order 1/m by the a matrices. It can be reduced to order 1/m 2 by 
another transformation 

ip — > exp(— iga ■ E/4m 2 )tp , etc.. (13) 
The non-relativistic lagrangian to order 1/m finally is 



£ N RQCD(x) = -\F* u Fr + ^l°D -m)^ + ^^ + ^^-^i) } (I IS 



which can be written in terms of uncoupled 2-component spinors by decomposing if> = 
(u, v*) i/i = (u\ —v) and using the Dirac representation of the 7-matrices. 

To effect a short distance expansion we separate the gluonic fields in slowly varying 
background- |TB[ (A) and quantum- (Q) fields having high fourier components: 



V v a = Al + Ql. (15) 
The division flTSP can however only be defined if the gauge is fixed. In deriving the 



spin-dependence of nonperturbative interactions for heavy quarkonia Curci et al. |L9] found 
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a particular gauge very convenient. We essentially follow their procedure here and impose 
the Coulomb background gauge condition 

AG* = o (16) 

for the quantum fields, where 

DpQ,, = d^Q v + gA^ xQ u = d^Q u + gf abc A^Q vc . (17) 



The background fields are defined in a modified Schwinger gauge JL9| 

A) = -\f*? ; K = -F^x\ (18) 

valid to order x 2 , where we assumed that the field-strengths corresponding to the background 
fields are constant (or have sufficiently low momenta, such that they can be regarded as 
essentially constant over the extent of the meson). This definition of the background fields in 
terms of (practically) constant field strengths also gives a definite meaning to the separation 
in equation (|15|). It also implies that the background fields in this gauge are (practically) 
time independent - a property which will become useful when a Hamiltonian is required. 

We next expand the nonrelativistic Lagrangian only to second order in the quantum 
fields and subsequently integrate them out in favor of an effective (coulombic) interaction. 
Although ghost terms are necessary in the gauge defined by (pJl), they do not contribute 
to quadratic order in the quantum fields. Our truncation of the interaction terms for the 
quantum fields eliminates all radiative corrections. To obtain them one would have to go 
beyond this approximation and calculate perturbative corrections before eliminating hard 
gluons. Fortunately, the asymptotic freedom of QCD guarantees that they are only logarith- 
mic at short distances and could be accounted for by a running coupling constant. These 
logarithmic corrections to the Coulomb potential do not seem to be dramatically important 
for describing heavy quarkonia spectra |3[ and we will not include them in this study. 

The matrix elements of the slowly varying background fields (A) will however be pa- 
rameterized. Their amplitude is large and an expansion in the coupling in this case not 
applicable. 

12 



After elimination of the quantum gluonic fields by their equations of motion the effective 
lagrangian in terms of the background and heavy quark fields becomes 



L eff = J d^xi—F^F^ + ^7°<9 + 91° A - m)ip 

2m 2m 
+g 2 J d 4 y^(x)T^(x)V ab (x, y)^{y)T b ^{y) (19) 

' r d 4 yJ?(x)(V? b (x,y) - K${x,y))j}(y) , 



where 

J a k (x) = -^-[(d k i>)Ty - $T a d k i> - ie ijk dST a ^)} 

+ ^(T b T a + T a T h )4>A a k - ^-f ahc e ijk A\^T c Y, 3 i, 

_ 2g 2fadc F d a J d 4 zV cb [x, z)^(z)T b ^(z) . (20) 

The propagator V relates the Q field to its source 

Q a (x) = Jd 4 yV ab (x,y)f b (y), (21) 



with 



f = g^T a ij + 2gf abc F m Q cl , (22) 



and satisfies 



{D j D j ) ab 1^ d {x, y) = 5 ad 5\x - y) . (23) 

Similarly the propagator T> appears in the elimination of the spatial components Qi and 
satisfies 

j d 4 zM« b (x, z)V%(z, y) = 5 ac 5\x - y)8 ik , (24) 
with the differential operator M? b given by 



13 



M$(x, z) = 5\x - z)[-(£> M Lg% - 2gf abc F* 

_|_^(T a T 6 + T h T a )i,5 l3 + ^f abc e ijk $T c Z k i>] (25) 
2m 2m 

+% 2 / d/a r e ^ / o(^)^ c (x, z)F^(z) . (26) 

Finally, the propagator /C enters when the lagrange multiplier of the gauge fixing condition 
for the quantum fields flTED is eliminated in turn. Its equation of motion is 

J d 4 zD^(x)V^(x, z)D dc {z)IC cd {z, y) = 5 bd 5{x - y) . (27) 

These rather formidable integro-differential equations for the Green functions can for- 
mally be solved order by order in the background field A. Since we will only retain terms 
of the hamiltonian matrix proportional to the lowest dimensional condensate < g 2 FF >, 



we only keep terms up to second order in the background fields in this gauge dig) . As we 
will see shortly, linear terms in the background field have to be retained, although we will 
assume that < F >= 0, e.g. that global colour- and lorentz- invariance is not broken. 

In order to obtain a tractable hamiltonian, further approximations are necessary. We will 
neglect all retardation effects in the effective interaction. This instantaneous approximation 
is correct to order 1/m, because retardation effects are generally expected to be of order 
1/m 2 . Since the interaction is instantaneous, the effective hamiltonian becomes time inde- 
pendent even in the presence of the background fields (which are (nearly) time independent 
in the gauge (fL8f)). This greatly simplifies the interpretation of our results. 

The effective hamiltonian, correct to order 1/m, r 2 and a s therefore is 



H = J d 3 x^u\x)mu(x) + v(x)mv\x) — -F^F^ 
-u\x)T A gEfx iU {x) - vfflfAgEfxiV^x) 



~* 2 —, 2 

-uUx)^u(x) - v(x)^vUx) 
v '2m y 1 w 2m K J 

+a J d 3 yu\x)T A u{x)\{y)f A v\y) (28) 
+ l^ k 9B C k f ACB J dhu\x)T A u{x) V -^v{y)T A v\y) 



f d 3 yu\x)T A u(x) { -^^v(y)T A vH 



64 J v ' r 
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1 

m 



7T 2 2 f 2 TT 2 (f) 2 f 2 

128 K ' K J 128 K J K ' 
■-u ] (x)T eij k gB k Xjd Xi u(x) - -v{x)T A e ijk gB k Xjd Xi v\x) 



-^(x^gBfuix) - \v(x)a{T A gB A v^ 



x) 



+ ™ f Anc gE n {dx ^ {S))T ^u{x)-u\x)T A {d x M{x)) 



-ie m d Xl {u\x)T A a k u{x) ■ / d 3 y(-36 tJ r + -±-L) v (#)T c V 



i a 



_^ f ADC gE D ^ dxtV ^))T A v \x)+v(x)T A {d Xi v\x)) 



J 



+ie ilk d Xl {v{x)T A a k v\x) ■ j d s y{-?>5 %J r + ^)u\y)T c u{y) 



Here u(x) and v(x) denote the anihilation operators for a quark and antiquark of mass m 
respectively whose spin and color indices have been suppressed, r = x — y and T A , T A are 
the hermitian generators of the SU(3) color Lie-algebra in the 3 and 3 representation. 



This is essentially the Hamiltonian to order 1/m derived previously by Curci et al. [19 



in position space, except for a term which can be regarded as a long range correction of the 
Coulomb potential 

3?rVa 



64 



d 3 xd 3 yv?(x)T A u(x) r 3 v(y)T A v 



(29) 



We have disregarded this term because it is of order r 3 and our gauge fixing condition ( 18|) 
is not valid at this order. 

The terms of the Hamiltonian ( p8|) linear in the chromo -electric and -magnetic fields as 
well as those proportional to x 2 are obviously also not translationally invariant and therefore 
apparently depend on the chosen gauge fixing point in fll8|) . This gauge dependence of (|28|) is 
however absent ]2D| in its matrix elements in the basis @,(§),(§D and (0). It cancels against 
terms which appear when derivatives act on the color transport matrix of the states. This 
cancelation of course only works to a certain order in the short distance expansion and only 
occurs if the Hamiltonian and the basis in which it is diagonalized are defined consistently. 
This is not a big surprise in any gauge theory, but a gauge invariant scheme to include 
nonperturbative aspects of the ground state in a hamiltonian formulation without losing 
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gauge invariance was only first proposed in |12], but not used to its full extent there. From 



the above we see that a short distance expansion in the construction of the basis as well as 
the hamiltonian to order r 2 can be performed and the effects from a nontrivial expectation 
value 4> 2 consistently included. 

The straightforward but lengthy computation of the matrix elements of the Hamiltonian 
(p8|) in the basis for pseudoscalar mesons (|) to obtain the coupled differential equations C|32|) 
for the wavefunctions ||2(J will not be exhibited here. We do however also have to account 



for matrix elements of the purely gluonic part of the Hamiltonian 

H G = J d 3 xH G (x), H G (x) = (E 2 + B 2 )/2 + 0(a) (30) 

in the basis states (0),(iD;(iD an d (0)- The order a terms in (|30[) are due to the 1-loop 
corrections from the quantum fields. When these corrections are properly included, the 
above Hamiltonian should be consistent with the trace anomaly |25j and reproduce the 



relation between the energy density of the nonperturbative gluonic ground state and the 
condensate value. All we will need in the sequel is that this correction is a local operator of 
dimension 4 and order a (such as (lla/167r)(i? 2 — B 2 )). 

From rotational symmetry we conclude that only diagonal matrix elements of H G can 
be non-vanishing. Since we neglect expectation values of operators with dimension greater 
than 4, matrix elements in "octet" -states (0) vanish in this approximation. 

From (§) we might naively infer that the matrix elements of Hq between "singlet" states 
vanish as well. Phenomenologically we do however need an energy splitting between the 



singlet"- and "octet"- states. Leutwyler [26 attributes it to an effective mass of the low 



frequency gluonic modes. A more careful examination of contributions from H G reveals 
the origin of such an effective mass in the present approach. To obtain the energy of the 
meson relative to that of the vacuum, one has to commute the purely gluonic hamiltonian 
with the creation operator for the meson. As discussed previously, gauge invariance requires 
that this creation operator contains a gluonic string, or color transport matrix. Within our 
approximation, we therefore have to consider the commutator of the creation operator for 
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the singlet component with the local Hamiltonian density of TCg(x). Defining the creation 
operator for the singlet component of the meson by 



the contribution of the lowest dimensional operator to the commutator is of the form 



where A(x) is a gauge invariant scalar operator of dimension 4 and order a (because the g°- 
part of the string obviously commutes) and 9y{x) is a dimensionless c- number function that 
vanishes for x outside the meson (because the string in Og and x are otherwise separated 
by a spacelike distance). Since we do not neglect the vacuum expectation of the gauge 
invariant operator aF 2 in our approach, the matrix element in the singulet channel within 
our approximation becomes 



where C is a constant which depends on the exact nature of the function 9y(x). 

Although we cannot determine the value of this constant on theoretical grounds, the 
above argument shows that it would be inconsistent with our approximations to neglect 
this contribution of Hq in the singlet channel. A non-vanishing constant C implies that the 
lowest order energies of the singulet and octet states differ by gluonic contributions. From 
the correlation length of (Q) discussed in the introduction one would expect a splitting of the 
octet and singulet channels of the order of 600 MeV, if the Fock-basis is a reasonable one. 
Our best fit to the quarkonium spectrum was obtained with C = —756. That this splitting 
of the "singlet" and "octet" Fock states is close to the correlation length of (|2]) indicates 
that the Fock-space is quite adequate for this calculation (i.e. the interaction energies are 
small compared to the splitting of the bare Fock states). 

m g = —C can also be interpreted as an effective mass for the low frequency gluon mode 
in the octet-channel (since a small change in the overall normalization of the energy can be 




[H G (x),Ot] = 0+A(x)9 v (x), 




(31) 
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absorbed in a slight change of the heavy quark mass). This is the point of view taken by 
Leutwyler ||26|| . Our estimate m q = — C ~ 756MeV is not in conflict with the fact that the 



lightest glueball has at least two such excitations and a mass of ~ 1700MeV ||10|| , although 
a somewhat smaller glueball mass would seem more natural. 

Including this energy shift between "singlet" and "octet" components, the coupled set 
of differential equations for flavor singlet pseudoscalar mesons becomes 



2m - E + C 



1 d 2 



7T— + 



TT 2 (b 2 r 2 



m dr 2 3 r 

, irdr q 
+ h= + 



36m J 5(r) 
g 2 (pr d 



+ 



3y/2 4 v / 2m %^/2m dr 



2m — E — 



1_ (?_ 

m dr 2 



2 . la s 
— ) H 



% 2 p / \ 
64m 2V ; 

3\/2 ' Sv^m 8V2mdr 



it<pr g 



)S(r) 



1 , 5 : 



1 a s 



2— + ~ 



6 r 



5^ 2 



64m 



5(r) 



(32) 



where E is the mass eigenvalue of the quarkonium and m is the mass of the constituent 
quarks. The functions S(r), E\(r) and -E^tj) are related to the wave function components in 
the expansion (||) via 

Mr) = $S(r) ; feW = MM ! = pr^(r). 

Note that there is no coupling to the pseudoscalar B-states, because it is proportional to 
the total spin in (p8|) , which vanishes for pseudoscalar mesons. As indicated above, these 
equations only depend on the relative coordinate r and all reference to a gauge fixing point- 
present in (p8|)-has disappeared in the evaluation of the effective Hamiltonian in the gauge 
invariant basis. 



Apart from the constituent quark mass, the only parameters that enter equations ( j32|) 
are the (reasonably well known) condensate <fi 2 and the constant C, whose value we have 
tried to estimate above, as well as the strong coupling constant a s . 



IV. "EFFECTIVE POTENTIAL" 



Before solving the coupled set of equations ( |3"2"D by numerical methods, it is instructive 
to consider the limit of infinitely heavy quarkonia. In this case, all terms proportional to 

18 



1/m (including the kinetic energy) in ( |32"D can be dropped, and the resulting equations give 
the binding energy V = E — 2m for states where the quark and anti-quark are localized 
a distance r apart (i.e. for wave-functions S(r), El(r) and E2(r) all proportional to 5{r — 
ro). It is natural to compare this mass-independent binding energy to phenomenological 
potentials and those extracted from the expectation values of Wilson loops in numerical 
lattice simulations. 

In this limit the i?2-component decouples and one has to solve the algebraic equations 

(-V + C-~)S{r) = -^E 1 (r) 

(-V + l^Mr) = ^S(r). (33) 

The eigenfunctions are obviously localized and the eigenvalue or effective potential, V(r), 
given by 



V(r) = J. ( l^ + jM-^. 2 ^fll + C^ (34) 
2 V 6 r V 4r 2 3 r 9 



Figure |T](a) shows this effective potential and its Coulomb part for parameters <fi 2 = 
(SQOMeV) 4 , C = —756 and a s = 0.39, which we found appropriate for charmonium. 

Although this potential is certainly no longer valid for r > 0.9fm, where the root in (|34"D 
becomes purely imaginary it does show a nearly linear behaviour for intermediate distances 
0.4fm < r < 0.8fm, with a correspondingly constant force of ~ 840MeV/fm, which compares 
favorably with a string tension of about 800 — lOOOMeV/fm extracted from recent lattice 
calculations of the potential p7j| . In Fig. |l](b) we compare our effective potential ( |34]) to that 
extracted from lattice data |27] and to the phenomenological potential used by ref |nj. It is 



perhaps also of some theoretical interest, that this potential closely resembles that derived 
from an "instantaneous" Gribov type gluon exchange [fj! 



VWO := -\f J i^^T/ 1 ' - , (35) 

in the limited range of interest r < .9fm with k = 500MeV. For comparison we also show 
the potential © in Fig 0(b). 
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It is encouraging that our rather crude approximations to the vacuum structure seem 
to qualitatively reproduce the potential for very heavy quarks at small distances. The 
analytical results of this section justify the numerical calculation of "octet" -components in 
heavy quarkonia which we now present. 



V. NUMERICAL RESULTS 

Going beyond the static approximation by including the kinetic energy of the quarks 
but still neglecting coupling terms of order 1/m in (|32|) , the quarkonium spectrum becomes 
discrete and the following coupled set of differential equations must be solved numerically 



2m- E + C 



Id 2 4 a 



S(r) = VEi(r) 



m dr 2 3 r J 3\/2 



1 / d 2 2 , la.i , , nor 



2m-E--(—--) + 

m or z r z or 



E 1 (r) = ^=S(r) . (36) 



This is essentially Leutwyler's |7j approximation to the problem, who obtained the pertur- 
bation of the Coulomb spectrum in 2 nd -order of and found that it is exceedingly large for 
canonical values of the condensate. Note that the Coulomb force is repulsive in the "octet" 
channel and 1/8-th in strength compared to the "singlet" channel (this is just the ratio of 
T®Tg in the two representations) and the instate still decouples in this approximation. In 
this approximation, vector- and pseudoscalar- quarkonia are furthermore still degenerate. 
This is expected, since the spin splitting is of order 1/m. It however is another nontrivial 
consistency check of our method, because the basis states for vector mesons are of course 
quite different. Nevertheless equivalent equations to (P^| ) result, if 1/m potential terms are 
neglected. 

The relative minus sign of the two off-diagonal coupling terms in ( |3~6"D shows that this 
Hamiltonian is not hermitian and that its eigenvalues will generally not be real. This effect 
is completely missed if the coupling is treated as a perturbation. To any finite order in 
perturbation theory the correction to the Coulomb spectrum is real. Perturbation theory 



does however show that a few (low lying) eigenvalues of (|3~6"D are real for a sufficiently small 
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0. For the canonical value 2 = (300 — 380MeV) 4 we find numerically 2 — 4 stable bound 
quarkonia in the pseudoscalar channel, depending on the heavy quark mass m. 

But even for these low states, the deviation of the exact correction to the coulomb 
eigenvalue from the 2 nd -order estimate is large for the charmonium and bottonium systems as 
shown in table |. We conclude with Leutwyler ]7J that perturbation theory is not applicable 
for canonical values of 0, but that the effective non-hermitian Hamiltonian equations ( |36|) 
do yield reasonable corrections in an exact solution. 

The value C = — 756MeV used by us was ajusted to reproduce the correct splitting 



between r\ c and r( c |28| when 1/m terms are included (see below). The corrections to the 

coulomb spectrum are however not extremely sensitive to the value of C once it is large 

enough. Using C = — 1400MeV instead would only reduce the nonperturbative contributions 

for the ground and first excited bottonium-states in table p] to 6 and 75MeV respectively. 

The fact that a perturbative evaluation in is not appropriate in (^) can also not be 

circumvented by including loop corrections (higher orders in a s than we have treated so 

far) to the perturbative coulomb potential. Titard and Yndurain |J{J recently proposed to 

modify the perturbative part of the interaction in the following manner 

Os_ _^ a 8 (l^) [1 + {ax + 7£/V2)q s (/x 2 )/7r] (3 a 2 (f^ 2 ) log r/i 
r r 2ir r ' 

where the appropriate constants for the SU(3) color group with 4 light quark flavors are 



Po = 8.33 and a\ = 1.47. The first term of (|37|) which contains a piece of one-loop 
radiative corrections was taken by Titard and Yndurain as an effective Coulomb potential 
and solved exactly. The effective coupling constant is defined as 

a s ^ 2 ) = [1 + ( ai + lE p /2)^l}a s (fi 2 ) . (38) 

7T 



The second term in (|37j ) was treated to first order in perturbation theory. A new scale 
parameter n was introduced which depends on the quarkonium system under consideration. 
Taking the effective Coulomb potential we obtain the deviation of the eigenvalues from 2 nd - 
order perturbation theory in shown in table [TI| for two sets of parameters in the bottonium 
system (still neglecting 1/m corrections). 
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The deviation of the exact correction from the perturbative estimate is reduced some- 
what, especially in the second case, but still far from negligible. In assessing the quality of a 
perturbative evaluation in this case, one should also keep in mind that the scale parameter \x 
of Ref. [[ftj was chosen in such a way that the 2 nd order correction in is precisely canceled 
by the correction terms to the effective Coulomb potential in (157). For the first set of pa- 



rameters this cancellation occurs for the ground state energy. The second set was chosen so 
that the splitting between the first excited state and the ground state is not affected to 2 nd 
order. This is obviously a quite arbitrary procedure which requires an additional parameter 
and furthermore does not cure the problem that perturbation theory simply does not apply 
(as table |I| clearly indicates). 

We therefore will not include these modifications to the Coulomb force in our discussion 
of the numerical solution to the full set of coupled equations fl32|). The inclusion of 1/m 
potential terms lifts the degeneracy of pseudoscalar- and vector- quarkonia and we restrict 
our discussion here to the pseudoscalar case. 

The instates now couple in, but generally have small (negative) norms, because the 
coupling is of order 1/m. We nevertheless solved the full set of coupled equations numerically, 
although a perturbative treatment of the .E^-state admixture would probably have been 
sufficient. The eigenvalues we obtained are summarized in table |T| and compared to those 



obtained by Eichten et al. |3T[ with the phenomenological funnel potential. Note that we 
only found 3 or 4 stable solutions in the charmonium and bottonium systems respectively. 
We used the same values for the quark masses and the coupling constant as Eichten et al. 



|31|.The constant C was chosen to reproduce the experimental splitting (not confirmed p8| ) 



between t] c and rj' c . It was not adjusted in the bottonium system. The gluon condensate value 



lc 

<>2 _ /"}«nA/IVan4 



(360MeV) we used is within QCD-sumrule estimates [29| for this nonperturbative 
quantity. All eigenvalues were finally shifted by E = 98MeV to give the correct 7] c mass. 
(This small shift in the overall energy normalization can be eliminated by a slight change 
of ~ 50MeV in the quark masses used by Eichten et al. |BT|] and a corresponding small 
adjustment of the other parameters. To have a more direct comparison of the wavefunctions 
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and spectra, we refrained from making these adjustments here.) 

In Figs. and |3] we show the eigenf unctions for the various components of our quarko- 



nium states (|9]) as well as the corresponding eigenfunction of Eichten et al. f3~I |. The singlet 
component of our ground state wave functions are very similar to those of the funnel po- 
tential. At small radii this is true also for the excited states, since the coupling to "octet" 
components is proportional to r in (^). The "octet" -components increase with increasing 
excitation energy of the quarkonium and lead to the appearance of extra nodes in the higher 
lying "singlet" wave functions at large radii (since this is a coupled channel problem, the 
extra nodes do not mean that we missed some bound states). As noted earlier, the E- 
components of the meson state contribute negatively to its norm. All the stable quarkonia 
states we found are however positive norm states. We could not obtain any stable state 
where the octet components are dominant. 

Let us speculate at this point on the fact that only very few stable quarkonium states 
were found. This is of course quite in line with the experimental observation that only a few 
heavy quarkonia are stable against decay by strong interactions. For reasons which we had 
not anticipated, this basic property seems already to be incorporated in the nonhermitian 
coupling to ^-components. 

The strength of this coupling in our model is however determined by the gluon condensate 

we see 



(j) , which is not expected to vanish even in the purely gluonic theory. From table |ITT 
that the instability sets in at an excitation energy of between 1 — 1.3GeV in this model. Since 
we cannot account for the decay into light gg-mesons with a parameter which is essentially 
independent of the number of (light) flavours, we speculate that the nonhermitian coupling 
proportional to effectively accounts for the possible decay channel 

Quarkonium* — > Quarkonium + Gluonium , (39) 

in this model. From the fact that we do not seem to find any stable quarkonium states 
more than 1.3GeV above the ground state, we would estimate this to be the threshold 
for the production of the lightest gluonium. This estimate of the lightest gluonium mass 
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me ^ 1.3GeV, is in almost perfect agreement with our previous interpretation of the energy 
shift — C = m g ~ uic/2. 

Since the production of light gg-pairs has a much lower threshold, it is this process which 
limits the stability of physical quarkonia. We therefore expect this model, which does not 
(not even effectively) incorporate this decay channel, to still predict more stable states than 
are experimentally observed. 



VI. CONCLUSION 

We developed a hamiltonian formalism, which enabled us to estimate the effect of a 
nonperturbative gluonic ground state on quarkonia in a systematic short distance, weak 



coupling and 1/m expansion of the effective hamiltonian. The gauge invariant basis [12] 
was extended to include color octet quark-antiquark pairs coupled to vacuum fluctuations. 
Hamiltonian matrix elements in this basis are gauge invariant to the order in the short 
distance expansion we considered. 



After separating hard and soft gluons in the gauge (0), we obtained the effective hamil- 
tonian neglecting radiative corrections to the coulomb interaction from hard gluons. We 
thus neglected the logarithmic corrections to the effective coupling strength at very short 
distances. The correct behaviour of the potential for r < 0.2/m can in principle be in- 
cluded by "hand" in a modification of the coulomb part of the interaction ||23|| . Although 
the quarkonium spectrum is not very sensitive to this correction at small distances, it could 
become important for the evaluation of decay widths, which depend on the wave function 
at the origin. 

In the limit of very heavy quark masses, where all terms of order 1/m can be neglected, 
an energy independent effective potential was obtained, which shows an approximately lin- 
ear behaviour at intermediate distances 0.4/m < r < 0.8 fm with an effective string tension 
of ~ 8A0MeV. In our approach this behaviour of the potential arises due to the nonpertur- 
bative structure of the gluonic vacuum parametrized by its gluon condensate and was not 
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assumed from the outset as in most phenomenological quark models. This potential com- 
pares favourably with recent lattice results pTfl , the discrepancies at very small r < 0.2 fm 
being due to our neglect of radiative corrections. Surprisingly, our potential is almost ex- 
actly reproduced by an instantaneous interaction derived from the effective gluon exchange 
proposed by Gribov |22|| . Our short distance expansion for the effective potential however 
is only valid for r < 0.9/m, beyond which the potential aquires an imaginary part. An 
extension of the model to larger distances would require a more detailed knowledge of the 
vacuum structure in the form of higher dimensional condensates, or some other effective 
parametrization of this structure. The approach in this case would become increasingly 
phenomenological and also more complicated in this case. Its predictive power is therefore 
probably limited to heavy quarkonia, where a detailed knowledge of the potential for very 
large radii does not seem necessary. 

We showed that the numerical solution of the coupled channel problem for vector- and 
pseudoscalar- quarkonia (they are degenerate to order 1/m) is feasible and an exact diago- 
nalization of the hamiltonian in the extended basis therefore possible. The resulting exact 
spectrum does not show the far too rapid increase of the eigenvalues with the principal quan- 
tum number of the perturbative approach to the vacuum effects |j] 0. In addition to the 
usual "singlet" wave-functions describing the quark and anti-quark of the quarkonium when 
they are coupled in a colour singlet, we also obtain the "octet" components of quarkonium 
states describing the quarks in the octet configuration when an additional (soft) gluon is 
around. Since hadronic decays mainly proceed from this "octet" configuration with the cre- 
ation of an additional octett gg-pair from a hard gluon, this approach opens the possibility 
of estimating nonperturbative contributions to hadronic decays. 

We compare our results for the spectrum and wavefunctions with those of the Cornell 
potential |3]J for pseudoscalar quarkonia in order 1/m. With the standard value for the gluon 
condensate and quark masses and coupling constant used by the Cornell group we obtain 
the correct splitting between t] c and r]' c and make predictions for the rjbS. Our main concern 
was however a better theoretical understanding and justification of the phenomenological 
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ingredients common to most nonrelativistic models for heavy quarkonia and we refrained 
from adjusting the few parameters of this approach to optimally reproduce the experimental 
data. A better description of the potential at short distances with the inclusion of radiative 
corrections to the coulomb force and the consideration of hadronic decay channels is clearly 
desireable before a detailed comparison with phenomenology is attempted. 
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FIGURES 

FIG. 1. (a) Effective potential (solid curve) from ( |34] ) and Coulomb potential (dotted curve) 
with a s = 0.39, C = -756 MeV and (j) 2 = (360 MeV) 4 . 

(b) Effective potential (solid curve) with the same parameters as in (a). The potential extracted 
from lattice data [^7]] with yfa = 365 MeV (dot-dashed curve) and yfa = 505 MeV (dashed curve) . 
Gribov potential |^] with k = 500 MeV (crosses) and Cornell potential |HJ (dotted curve) with 
a s = 0.39 and a = 2.34 GeV" 1 . 

FIG. 2. The wavefunctions of the ground, 1 st - and 2 nd - excited pseudoscalar states of char- 
monium are shown in figures (a), (b) and (c) respectively. The dashed curve is the wavefunction 



for the funnel potential [31] for comparison. The solid curve is the singlet component S(r) of the 
quarkonium state in our calculation. The dot-dashed and dotted curves are the "octet" compo- 
nents E\(r) and i?2( r ) of equation ( j32|) respectively. The solution was obtained with the parameters 
m c = 1840 MeV, a s = 0.39, <p 2 = (360 MeV) 4 and C = -756 MeV. 

FIG. 3. The wavefunctions of the ground, 1 st -, 2 nd - and 3 nd excited pseudoscalar states 
of bottonium are shown in figures (a), (b), (c) and (d), respectively. The dashed curve is the 
wavefunction for the funnel potential |?l| for comparison. The solid curve is the singlet component 
S(r) of the quarkonium state in our calculation. The dot-dashed and dotted curves are the "octet" 
components E\(r) and E2(r) of equation (^) respectively. The solution was obtained with the 
parameters m b = 5170 MeV, a s = 0.39, <p 2 = (360 MeV) 4 and C = -756 MeV. 
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TABLES 

TABLE I. The Coulomb energy in MeV is presented in the first column for the ground state 
and first excitation of cc and bb (pseudoscalar or vector). Second and third columns contain the 
nonperturbative contributions in MeV calculated respectively within perturbation theory and with 
our method. We used m c = 1840 MeV, m b = 5170 MeV, a s = 0.39, (j) 2 = 0.012 GeV 4 and 

C = -756 MeV. 

Coulomb Pert, theory Non-pert. 

ric, J/Y> -124.4 311 66 

rj c , </>' -31.1 21525 424 

rjb, T -349.5 14 8 

rf b , V -87.4 970 111 



TABLE II. The Coulomb energy in MeV is presented in the first column for the ground state 
and first excitation of bb (pseudoscalar or vector) at two different scales /x. Second and third columns 
contain the nonperturbative contributions in MeV calculated respectively within perturbation the- 
ory and with our method. We used C = -756 MeV, <f> 2 = 0.042 GeV 4 and for \i = 1.44 GeV: 
m b = 4866 MeV, a s = 0.38. For p = 0.99 GeV: m b = 5010 MeV, a s = 0.54. 



Coulomb Pert, theory Non-pert. 
rj b , T (fi = 1.44 GeV) -312 25 11 
rf b , V {n = 1.44 GeV) -78 1762 129 
rjb, T (n = 0.99 GeV) -649 6 3 
rf b , V {ii = 0.99 GeV) -162 396 72 
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TABLE III. Masses in MeV of pseudoscalar quarkonia with the funnel potential and with 



our effective nonperturbative hamiltonian. Parameters used: m c = 1840 MeV, ra b = 5170 MeV, 
a s = 0.39, <j) 2 = (360 MeV) 4 , C = -756 MeV. 

funnel nonpert. 

f] c 2980 2980 

rf c 3571 3594 

r# 3994 3993 

7] b 9213 9344 

rf b 9805 9739 

7]'l 10150 10084 

n'l' 10427 10610 
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